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Abstract 

We present a time-dependent localized Hartree-Fock density-functional linear response approach for the 
treatment of photoionization of atomic systems. This approach employs a spin-dependent localized Hartree- 
Fock (SLHF) exchange potential to calculate electron orbitals and kernel functions, and thus can be used to 
study the photoionization from atomic excited states. We have applied the approach to the calculation of 
photoionization cross sections of Ne ground state. The results are in agreement with available experimental 
data and have comparable accuracies with other ab initio theoretical results. We have also extended the 
approach to explore the photoionization from Ne excited states and obtained some new results for the 
photoionization from outer-shell and inner-shell excited states. 
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I. INTRODUCTION 



Making use of local potentials and independent-particle response functions, density functional 
theory (DFT) combined with linear response approximation (LRA) 1, 3] has been successfully 



applied to study dynamic processes, such as photoabsorption P, 



6, 



polarizability 



11 



and dynamic 



13l.ll4l. Il5l|. of atomic and molecular systems. The most attractive features of 



such an approach are its satisfactory accuracy and computational simplicity and efficiency. How- 
ever, since the conventional DFT using traditional exchange-correlation (XC) potentials obtained 
from uniform electron gas, such as local density approximation (LDA) [igI. fl?! and generalized gra- 
dient approximation (GGA) 17. is a ground-state approach, the conventional DFT-LRA 
approach can only be used to investigate the dynamic processes associated with the ground state 
of a system. Even so, the counterpart of the DFT-LRA approach for the excited states has not yet 
been reported. 

The difficulty encountered in the extension of the conventional DFT-LRA approach to the 
excited states stems from the XC potential used to characterize the excited states. A qualified 
XC potential for the excited states is required to be symmetry-dependent and self-interaction free 
and have a correct long-range behavior. The symmetry (such as electronic configuration, electron 
orbital angular momentum, and electron spin) of the XC potential is used to distinguish a state from 
the others, the self-interaction-free property of the XC potential is used to make the calculation of 
electron orbital energy accurate, and the correct asymptotic behavior of the XC potential is used 
to guarantee the Rydberg virtual orbitals which play a key role in autoionization resonances Q]. 

Recently, a localized Hartree-Fock (LHF) density-functional approach has been proposed and 
successfully applied to the ground-state calculation of atomic and molecular systems j2ci]. In this 
approach, the LHF exchange potential is self-interaction free and exhibits the correct long-range 
behavior. It only needs occupied orbitals and depends on the orbital symmetry of the state. More 
recently, a spin-dependent localized Hartree-Fock (SLHF) density-functional approach has been 



developed for the excited-state calculation of atomic and molecular systems 
together with Slater's diagonal sum rule 



21| . This approach 



22| | have been successfully used to calculate the energies 



of multiply excited states of valence electrons of atomic systems 



excited states of closed-shell [23] and open-shell [2j| atomic systems. 



21[ and the energies of inner-shell 



In this paper, we present a time-dependent localized Hartree-Fock density-functional linear 
response approach for the treatment of photoionization from atomic excited states by combining the 
SLHF DFT with LRA. In this approach, the SLHF exchange potential is employed to calculate both 



2 



the Kohn-Sham (KS) electron orbitals and kernel functions. This is different from the approach 

n 

used in Ref. [10] where the LHF potential was only used to calculate the electron orbitals. We 
have applied this approach to the calculation of photoionization cross sections (PICS) of Ne ground 
state and extended it to the computation of PICS of Ne excited states. 



II. THEORETICAL METHODOLOGY 

A. Linear response approximation of photoionization 

Suppose that an atomic system is in a time-dependent external field along z axis E (t) = Sqc'^^^, 
where £q is the amplitude and uj the frequency. In dipole approximation, the interaction potential of 
an electron and the external field is (/>^^* i^,t) = <P^^^ ('^j'^) e"*"^*, where (j)^^^ (r, w) = z£q. Presence 
of the external field will induce a perturbation to the system and produce a redistribution of 
^. In the case of weak field considered here, the dynamic response of the system 



electron density 



has the same time dependence e as the field [8:] and thus can be described by the LRA 25l |. 



In frequency domain, total perturbingpotential that an electron experiences can be expressed, in 
spin-dependent DFT framework 3, g, 7, ^, as 

</>f (r, u:) = <A^-* (r, ^) + > J / K,.> (r, r') 5p,, (r', ^) dv\ (1) 



where, the second term on the right-hand side (RHS) is the field-induced potential, a is the electron 
spin (cr = t for spin-up or cr = | for spin-down), bp^ (r, w) is the field-induced electron density 
which is the deviation of the perturbed electron density from the unperturbed electron density 
Pa (r), and -fCcro-' (r, r') is the kernel function 

K,^, r, r' = + . 2 

^ ' |r - r'l bpa' (r') 

Here, the first and second terms on the RHS represent the field-induced changes of Hartree potential 

and XC potential Vxca (r), respectively. 

The field-induced electron density bp^j (r,a;) is related to the total perturbing potential (\)^'^^ 

by 

8p, (r, oj) = Jxa (r, r', (r', u;) dr' , (3) 



where, Xa is the complex susceptibility given by 

occ all 



J 



huj - {eja- - Eia) + ie hw + (sja - Eia) + ie 



(4) 



Here, 7^0- (r, r') = ip*^ (r) ipi^j (r'), ipi^ (r) is the ith. electron spin-orbital, e is a positive infinitesimal, 
and the notations all and occ represent that the sums run over all and occupied electron orbitals, 
respectively. 

From Eqs. ([1]) and ([3]), <j)^^^ (r, uj) and dp^ (r, uj) have to be calculated in a self-consistent field 
procedure. Applying Eq. (P) to Eq. ^ and introducing a kernel function 

N^^i (r, r', io) = j Xcr (r, r", u) K^^i (r", r') dr" , (5) 

one obtains an equation for 5pcr (r, uj) 

y [^cra'^ (r - r') - N„„: (r, r', w)] dp^i (r', dr' = j Xa (r, r', (^""^'^ (r, w) dr. (6) 

cr' 

The polarizability a (cu) is the ratio of the induced dipole moment to the external field strength 
a{uj) = -^X] j z5pa{Y,uj)dr 



The PICS a {uo) is calculated by 

a [lo) = Im a (w) , (8) 

where, c is the speed of light. 

The LRA is usually referred to as a time-dependent method since it takes the time-dependent 
field- induced electron density into account. In contrast, if the field-induced electron density is 
neglected in Eq. ([T|), one has (j)^^^ = 4>^a^- I^i this case, the cross section Eq. ([7]) is reduced to the 
result of independent-particle approximation, which is referred to as a time-independent method. 

B. Electron spin-orbitals 

?o calculate the occupied electron spin-orbitals we use the SLHF density-functional approach 
■'■^ this approach, the electron spin-orbital (^jg- (r) and orbital energy are calculated 



21 



from the KS equation 

r 1 _ „ 1 

ipia (r) = eiaipia (r) , (9) 



where, is the local effective potential given by 

^.^'^(r) = -7 + E/^^r' + y_(r), (10) 



In Eq. ([To]) , the spin-dependent electron density is calculated by (r) = Yld=i \ '^i<T (r)P, where 
No- is the number of electron with spin a. On the RHS of Eq. (jlOp . the first term is Coulomb 
potential, the second term is Hartree potential, and the third term Vxca- (r) is the XC potential. 



The XC potential can be decomposed into exchange potentia. 
Vctj (r). The SLHF exchange potential V^l'^^{r) is given by 

1 



21 



(r) and correlation potential 



23|, 



24 



26| 



V 



SLHF 



(r) 



+ 



P<T(r) 

-^|7.(r,r')l4r"(r')7.(r',r)rfr' 
la (r,r')7f7 (r',r")7<7 (r",r 



(11) 



where, 7o- (r,r') = Yld=ili<r {'^i'^')- the RHS of Eq. (fTTj) . the first term is Slater potential [22| 



and the second and third terms are the corrections to the Slater potential. The V^^J'^^(r) behaves 
asymptotically as Slater potential and thus approaches to the correct — 1/r at long range 26l |. 
As for the correlation effect, the second-order gradient correlation potential and energy functional 
proposed by Lee, Yang, and Parr (LYP) [27] can provide an excellent correlation energy for atomic 
systems and will be incorporated into the calculation to estimate the correlation effect in this work. 



C. Green function 



In Eq. dH) the sum over j needs all the (occupied and unoccupied bound and continuum) 
electron orbitals. This makes it extremely difficult to accurately calculate X(T(r,r',u;) directly 
from Eq. dH). To circumvent this difficulty, a Green function associated with the KS equation has 

nnn 

been introduced to calculate [SjlajlSl- The Green function {r,r',E) is calculated by 



E+lv'-Vf{r)^ {r,r',E)=6{r-r'), 



(12) 



under appropriate boundary conditions. The Green function can be expanded in terms of a com- 
plete set of KS electron spin-orbitals as [s] 



-y* (r r') 



(13) 



Applying Eq. ^3^) to Eq. 1^ one has 

occ 

Xa (r, r', cj) = ^ [7,^ (r, r') (r, r', e^, + fuv) + 7*, (r, r') G^ (r, r', - hw)] . (14) 

i 

Thus with assistance of the Green function, only the occupied orbitals are needed to compute Xa- 



III. COMPUTATIONAL DETAILS 



A. Electron spin-orbitals 



The electron spin-orbitals of an atomic system can be calculated by using the procedure previ- 
ously developed in [21I, 2^]. In spherical coordinates, the electron spin-orbital ipi^ (r) is expressed 



as a product of a radial spin-orbital Rniaif) and a spherical harmonic Yim{6, 

Rrdu{r) 



•^ia (r) 



-Yy 



Im \ 



(15) 



where, n is the principal quantum number, I is the orbital angular momentum quantum number, 
m is the azimuthal quantum number, and i is a set of quantum numbers except the spin a. The 
radial spin-orbital Rnia{r) is calculated from the radial KS equation {21!. [2^ 



R 



nla 



(16) 



2 dr2 2r2 

where, v^{r) is the radial effective potential [21]. To obtain high-precision electron spin-orbital 
and orbital energy, we use generalized pseudospectral (GPS) method 2g] to discretize the radial KS 
equation (jl6p . The GPS method associated with an appropriate mapping technique can overcome 
difficulties due to singularity at r = and long-tail at large r of the Coulomb interaction. It allows 
for nonuniform and optimal spatial discretization with the use of only a modest number of grid 
points. It has been shown that the GPS method is a very effective and efficient numerical algorithm 



for the high-precision solution of KS equation 



21 



, I29I 



130(. 



B. Green function 



In spherical coordinates, the Green function can be expanded in terms of partial waves as 

(r, r', E)=Y^ {6, ^) Qlo (r, r', E) Ylm {9', , (17) 



LM 



where, Qlct {r, r', E) is the radial Green function, which, with the appropriate boundary conditions, 
is determined by an inhomogeneous equation 



2 or or 2 r^ 



GLa {r,r',E) 



5 {r — r') 



(18) 



Alternatively, the radial Green function can also be constructed from the solutions of a homo- 
geneous equation 

£ L{L + 1 



dr"^ 



+ 



+ 2vf (r) - 



ULka [r] = 0, 



(19) 
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where, k = \f2E. If ^Lfco- ij) is the solution of Eq. (fT9l) being regular at the origin and V'Lfco- (j^) 
the solution behaving asymptotically as rh^^\kr) (where h^^^ is the spherical Hankel function of 



m 



the first kind), the radial Green function can be calculated by 

g..(r,/,i^) = -| ^^^-(^<)t"'-(^>\ (20) 

where, r< (r>) refers to the smaller (larger) of r and r' and W = (pLkaip^ka ~ ^'ika'^Lka is the 
Wronskian oi (pLka and V'Lfccr- 



C. Absorber 

In principle, the boundary conditions at r ^ oo is required to calculate ipLka (?")• In reality, the 
boundary can not be set at r — > oo. No matter how far the boundary is, as long as it is located 
at finite distance, the out-going wavefunction Vl/co- {f) with £^ > may reflect on the boundary, 
making the PICS oscillating artificially. To remove the refiection we introduce an absorber for 
each out-going wavefunction. The absorber is characterized by an absorptive potential with a 
linear dependence of the radial coordinate 

0, < r < Ta 

jj r -Ta . ^ ' (21) 

-Uq — , ra<r< rmax 



U{r) 



where, rmax is the radial coordinate of the boundary, Uq and are two parameters representing 
the strength and starting position of the absorber, respectively. Obviously, Tmax — represents the 
width of the absorber. Similar absorbers have been used in the wavepacket method of molecular 
collisions 



31 



32( 1 and photoionization of molecules and atomic clusters recently [33|. It is shown 
that the PICS are insensitive to the absorber parameters. 

When taking the absorptive potential into account the out-going wavefunction is calculated 
from an equation obtained by replacing (r) with v'^ (r) + iU (r) in Eq. (jl9p . Since the behavior 
of an out-going wavefunction depends on electron spin-orbitals through both v'^ (r) and k{E) the 
absorber parameters may be different for different electron spin-orbitals. 

D. Susceptibility and cross sections 

In spherical coordinates, the susceptibility Xa (r, r', lo) can also be expanded in the partial waves 
HI 



X„ (r, r', u:)=Y, (0, <t>) Xla (r, /, co) Yi^ {9', . (22) 
Im 



From Eqs. ([H]), p5]) . (fT7|) . and (j22]) . the partial wave susceptibility x/o- {r-,r' ,uj) is calculated by 
Xl.{r,r,u) = ^2^2^Wn'i'a \{mO\LO)\ 

n'V L 

[Qlu {r, r', En'l'a + ^) + G*La {r, r' , En'l'a " ^)] • (23) 



Furthermore, expanding 5pa (r, uj) in the partial waves 6pa (r, uj) = X^Zm ^Pimo- (r, w) ^/m. 
using z = ^J^rYiQ {0,(j3), and from Eq. ([6]), we obtain 



+ ^ / Ni^a' {r,r',uj) 5pima' {r',u}) r''^dr', 



(24) 



where, Ni^^^' {r,r',uj) is the partial wave component in the partial wave expansion iVo-o-' (r, r') = 
Ylim ^im ^) -^io-cr' {r, r', uj) Yim {0' , (j)')- Solving Eq. (f24|) one obtains the partial wave component 
5 Pima (^i"^)- From Eq. ([7]), the polarizability a {to) is given by 

a (w) = - ^ V / (5pioa (r, cj) r^dr. (25) 
to V 

The PICS a (lj) is calculated from Eq. ([8]) by using a (w). 



IV. RESULTS AND DISCUSSION 

A. Photoionization from the ground state of Ne 

To test the approach developed in the preceeding sections we first apply it to the calculation of 
PICS from the ground state of Ne. The total PICS is plotted in FIG. 1, where the solid and dashed 
lines are the results of time-dependent spin-dependent localized Hartree-Fock (TDSLHF) method 
and time-independent spin-dependent localized Hartree-Fock (TISLHF) method, respectively. Also 
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35l | for comparison. It shows that the 



plotted in this figure are recent experimental data 
TISLHF PICS are substantially larger than the experimental results in a broad range of lower 
photon energy. While the TDSLHF cross sections are much closer to the experimental results and 
show significant improvement in spite of a little bit underestimate near the ionization threshold 
(IT) of 2p electron. In the range of photon energy shown in FIG. 1, apart from a small stepwise 
enhancement due to the contribution of 2s electron photoionization at the IT energy of 2s electron 
46.405 eV, the TISLHF PICS are structureless and decrease monotonically above the IT of 2s 
electron. The TDSLHF PICS however appear a series of sharp resonances in the energy region 



8 



from 43.0 to 46.4 eV. In the higher energy region, both TDSLHF and TISLHF results are getting 
closer and merge to the experimental results. 

The sharp resonances in the PICS are produced by the resonant transitions 2s np from the 
ground state ls^2s^2p^ to autoionization states ls^2s2p^np with n = 3 to 10. In FIG. 2(a) we 
show the details of the resonances in open circles. To quantitatively analyze_the resonances and 
estimate the autoionization widths we fit the resonances to the Fano profile 



a {uj) = fJo 



r] 



1 + k2 



(26) 



where, k = 2 (w — Er) /F, Er is the resonance position, q is the profile index, F is the line-width 
of the resonance profile, rj'^ is the correlation coefficient, and uq is the cross section without corre- 
lation. In FIG. 2(a) we show the fitted resonances in the solid line. The fitted resonance profile 
parameters are given in TABLE [J in the rows with TDSLHF. Also listed in TABLE [J are the avail- 
able experimental results [37] and the theoretical results of the time-dependent density functional 
method employing exact exchange-only KS potential and adiabatic-local-density-approximation 
exchange-correlation kernel with shifted 2s orbital energy (S-EXX/ALDA) 10 , time-dependent 
local spin-density approximation with shifted 2s orbital energy (S-TDLSDA) [6 1, time-dependent 
local density approximation (TDLDA) Q], and R-matrix method (R-matrix) [38], for comparison. 
It is shown that the line-width F and cross section ctq of TDSLHF are in overall good agreement 
with the experimental values and better than the S-EXX/ALDA, TDLDA and R-matrix results. 
The line profile index q of TDSLHF is larger than both the experimental and theoretical results. 
Except the S-EXX/ALDA and TDLDA results, the correlation coefficient rj'^ of TDSLHF is smaller 
than the others. As for the resonance position Er, all the resonances in TDSLHF PICS are about 
2 eV shift to the lower photon energy with respect to the experimental data. Similar phenomena 
were also found in P, 1^. This discrepancy mainly stems from the deviation of the 2s electron 
orbital energy from the experimental value [6|]. For demonstration, we notice that the 2s electron 
orbital energy of TDSLHF is —1.707 a.u. which is about 2.04 eV higher than the experimental 
value —1.782 a.u. [s^. To explore influence of the orbital energy to the PICS and resonance 
profiles we have also performed a TDSLHF calculation with the 2s electron orbital energy being 
replaced by the experimental value in the calculation of the Green function. The calculated and 
fitted PICS in the resonance energy region are shown in FIG. 2(b) in the open circles and solid 
line, respectively. The fitted resonance profile parameters are listed in TABLE |T] in the rows with 
S-TDSLHF. It is shown that with the 2s electron orbital energy replaced by the experimental value 
the resonances shift to the higher photon energy about 2 eV and the calculated resonance positions 



Er are in very good agreement with the experimental results. In addition, all the other resonance 
profile parameters, particularly F and rj'^, are also significantly improved and much closer to the 
experimental values. 

B. Photoionization from outer-shell excited states of Ne 

As an extension we apply the proposed approach to the computation of PICS of Ne outer- 
shell excited states. In FIG. 3, wc show, in the solid and dashed lines, the total PICS of TD- 
SLHF method and TISLHF method for the photoionization from the outer-shell excited state 
lsj^ls-f2sj^2s|2p|2p|3s| of Ne, respectively. In FIG. 4 we show the details of autoionization reso- 
nances in the total PICS. It is shown that the total PICS can be divided into five regions A~E. 
They are separated by the steps at ITs in the TISLHF PICS. The first region A starts from the 
IT of 3s I electron and ends at the IT of 2p^ electron followed by two series of sharp autoionization 
resonances in the TDSLHF PICS. These two series of resonances are produced by the transitions 
of a 2j9| electron to the higher bound S| and d-^ orbitals, respectively, and thus converge to the IT 
of 2p| electron, as shown in FIG. 4(a). Unlike the PICS of the ground state which increase with 
photon energy near the IT of 2p electron, the PICS of the excited state decrease with photon energy 
near the IT of the 3s-f electron. The second region B starts from the IT of 2p| electron and ends at 
the IT of 2pi electron. In this region there are also two series of sharp resonances in the PICS, as 
shown in FIG. 4(b). They are induced by the transitions of a 2p| electron to the higher bound 
and di orbitals, respectively, and converge to the IT of 2p^ electron. The third region C starts from 
the IT of 2pj electron and ends at the IT of 2s^ electron followed by a series of sharp resonances 
produced by the transitions of a 2s| electron to the higher bound p| orbitals, as shown in FIG. 
4(c). Note that one resonance produced by the transition 2s| 2p| in this series is located in the 
lower energy region A in FIG. 4(a). The fourth region D starts from the IT of 2s| electron and 
ends at the IT of 2si electron. This region consists of a series of sharp resonances produced by the 
transitions of a 2s ^ electron to the higher bound p^ orbitals as shown in FIG. 4(d). The last region 
E covers the energy region from the IT of 2sj electron to the maximum energy shown in FIG. 3. 
In this region the PICS are structureless since the next autoionization resonances occur at very 
high energy region when a Is-f electron is resonantly pumped to the higher bound p-^ orbitals. We 
have also performed the calculation of PICS from other outer-shell excited states of Ne. The total 
PICS have the similar structures as those from the excited state lsj^ls|2sj^2s|2p|2p^3s| shown in 
FIG. 3. 
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For the photoionization from Ne excited states, particularly for those with autoionization res- 
onances, both experimental and theoretical results are scarce. Thus it is hard to make direct 
and comprehensive comparison of our results to the experimental and other theoretical results. 
Through a roughly qualitative comparison we found that for the photoionization from the excited 
state lsjls|2sj2s|2p'^2p^3p^ the results of TDSLHF method are in the same order of magnitude 
as the available experimental data and other theoretical results 



41 



42( 1 . For example, at 



the photon energy 3.41 eV, the PICS of TISLHF and TDSLHF methods are 3.99 Mb and 6.70 
Mb, respectively, for the photoionization from the excited state ls|ls|2s|2s|2p|2p|3p|, while the 
experimental results for the photoionization from the excited state ls^2s^2p^3p "^D-^ is 2.15 Mb 



-field approximation for the photoionization from the ex- 



40| , the theoretical results of the centra 
cited state ls^2s^2p^3p is about 5.0 Mb [4l|, and the theoretical results of the single-configuration 
Hartree-Fock method for the photoionization from the excited state ls^2s^2p^3p is about 4.5 
42 1 . For the photoionization from ls^ls^2s i2si^ 2p'^2p'^3si , the results of TDSLHF method are 



Mb 



41 



43i,y,y. 



larger than the available experimental [431] and other theoretical results 

Interchannel interference has a significant impact on the autoionization resonances in the pho- 
toionization from excited states. Due to the interference of autoionization channels, the profiles of 
autoionization resonances, particularly those produced by the transitions from the same electron 
orbital such as those shown in FIG. 4(a) and (b), are no longer the standard Fano profile. Thus 
it is impossible to accurately evaluate autoionization resonance profile parameters by fitting the 
resonances to Fano profile. However, we can still estimate the peak position for each autoioniza- 
tion resonance very well. In TABLE [III we list the peak positions of autoionization resonances for 
the photoionization from Ne outer-shell excited states Isils-i2si2si2p'^2p'^nli^ with n = 3 ~ 4 and 
/ = ~ 2, where {} = ls^ls|2sj2s|2p^2p| is the abbreviation of the core configuration. 



4i- 



One of the important physical processes during the photoionization is orbital relaxation 
Due to the orbital relaxation the peak position of an autoionization resonance is different from the 
orbital energy difference between the two unperturbed electron orbitals involved. For example, for 
the excited state l,Sjls|2sj2s| 2p|2p|3s|, the orbital energy difference between the unperturbed 
2p| and 2s-^ electron orbitals is 23.484 eV. While the peak position of the autoionization resonance 
produced by the transition 2s-|- 2p^ is at 22.994 eV, which is about 0.490 eV shift toward the 
lower energy with respect to the unperturbed orbital energy difference. Our calculation shows that 
the relaxation effect is larger for the autoionization resonance produced by the transition from a 
deeper inner-shell electron orbital to a lower-lying autoionization electron orbital. 
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C. Photoionization from inner-shell excited states of Ne 

As another application we extend the proposed approach to the calculation of PICS of Ne 
inner-shell excited states. In FIG. 5 we plot the total PICS of TDSLHF and TISLHF methods in 
the solid and dashed lines, respectively, for the photoionization from Ne inner-shell excited state 
Is'ilsi2s^2pppl3si. The PICS can be divided into four regions A~D. Each region starts from the 
IT of an electron orbital and ends at the IT of next lower-energy electron orbital. It is shown that 
the TISLHF PICS are structureless apart from a stepwise enhancement at each IT. The TDSLHF 
PICS, however, contain several series sharp autoionization resonances, as shown in FIG. 6. 

The region A starts from the IT of 3sj^ electron and ends at the IT of 2p^ electron followed 
by two series of autoionization resonances, as shown in FIG. 6(a). The two series resonances are 
produced by the transitions of a 2pi electron to the higher bound S| and di orbitals and converge 
to the IT of 2pi electron. The cross sections near the IT of 3sj electron again decrease with photon 
energy. The region B starts from the IT of 2pj electron and ends at the IT of 2p^ electron. In this 
region there are also two series of autoionization resonances. They are produced by the transitions 
of a electron to the higher bound Sf and orbitals and converge to the IT of 2p-^ electron, as 
shown in FIG. 6(b). The region C begins with the IT of 2p| electron and ends at the IT of 2s-f 
electron. In this region, there is a series of autoionization resonances induced by the transitions of 
the 2s| electron to the higher bound p^ orbitals, as shown in FIG. 6(c). The last region D starts 
from the IT of 2s| electron and ends at the maximum energy shown in FIG. 5. In this region, the 
PICS are structureless since the autoionization resonances will occur at very high energy. 

We have also estimated the peak positions for the autoionization resonances. In TABLE 
we list the peak positions of autoionization resonances for the photoionization from Ne inner-shell 
excited states ls|lsj2s|2p'|2p|n/| for n = 3 ~ 4 and Z = ~ 2. Our calculation shows again that 
the relaxation effect is larger for the autoionization resonance induced by the transition from a 
deeper inner-shell electron orbital to a lower-lying autoionization electron orbital. 

V. CONCLUSION 

In this paper, we propose a time-dependent localized Hartree-Fock density-functional linear 
response approach for the treatment of photoionization of atomic systems. In this approach, the 
occupied electron orbitals are calculated by solving the KS equation with SLHF exchange potential, 
the complex susceptibility is calculated by using the occupied electron orbitals and corresponding 
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Green functions, and the PICS are calculated by using the susceptibility. The relaxation of electron 
orbitals is taken into account through the Green functions. To remove the possible reflection of 
the wavefunction on the boundary an absorber is imposed for each out-going wavefunction in the 
calculation. The absorber is characterized by an absorptive potential with a linear dependence 
of coordinates. Since the SLHF exchange potential has a correct long-range behavior and can be 
used to accurately calculate the electron orbitals of atomic excited states, the proposed approach 
is suitable for the calculation of photoionization from excited states of atomic systems. We have 
applied this approach to the calculation of PICS of Ne ground state. The results are in agreement 
with available experimental and have comparable accuracies with other ab initio theoretical results. 
We have also extended this approach to the computation of PICS of Ne excited states. The total 
PICS can be divided into several regions. Each region starts from the IT of an electron orbital and 
ends at the IT of the next lower-energy electron orbital. In each region there is one or two series 
of autoionization resonances produced by the resonant transitions of an electron from inner-shell 
electron orbitals to autoionization electron orbitals with higher energies. The orbital relaxation 
effect is larger for the autoionization resonance produced by the transitions from deeper inner-shell 
electron orbitals to lower-lying autoionization electron orbitals. 
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Figure Captions 

FIG. 1 (Color online) Total photoionization cross sections from tlie ground state of Ne. Tiie 
solid and dashed lines are the results of TDSLHF and TISLHF methods, respectively. The open 
and solid circles are the experimental data of Chan et al. [34] and Samson et al. [s^. 

FIG. 2 (Color online) Autoionization resonances in the total photoionization cross sections near 
the ionization threshold of 2p electron, (a) The results of TDSLHF method and (b) the results 
of TDSLHF method with the 2s electron orbital energy being replaced by the experimental value. 
The resonance peaks 1 to 8 are produced by the transitions 2s np with n = 3 to 10, respectively. 
The open circles are the numerical results and the solid lines are the fitted results to the Fano 
profile. 

FIG. 3 (Color online) Total photoionization cross sections of Ne outer-shell excited state 
ls|ls|2sj2s|2j?'|2p^3s|. The solid and dashed lines are the results of TDSLHF and TISLHF meth- 
ods, respectively. 

FIG. 4 (Color online) Autoionization resonances in the total photoionization cross sections of Ne 
outer-shell excited state ls|ls|2s|2s|2p|2p|3s|. (a) The resonances produced by the transitions 
2p| ns| and n'd| with n > 4 and n' > 3. Note that the resonance produced by the transition 
2s| 2p| is embodied in this region, (b) The resonances produced by the transitions 2p^ ns| 
and ndj with n > 3. (c) The resonances produced by the transitions 2s^ — > np-^ with n > 2. Note 
that the resonance for n = 2 is located in the energy region in (a), (d) The resonances produced 
by the transitions 2s| — > npi with n > 3. 

FIG. 5 (Color online) Total photoionization cross sections of Ne inner-shell excited state 
ls|ls|2s|2p|2p|3s|. The solid and dashed lines are the results of TDSLHF and TISLHF methods, 
respectively. 

FIG. 6 (Color online) Autoionization resonances in the total photoionization cross sections of 
Ne inner-shell excited state ls|lsj^2s|2p|2p|3sj. (a) The resonances produced by the transitions 
2pi — > ns| and n'di with n > 4 and n' > 3. (b) The resonances produced by the transitions 
2pi ns| and ndf with n > 3. (c) The resonances produced by the transitions 2s| np^ with 
n > 3. 
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TABLE I: Profile parameters of the autoionization resonances produced by the transitions from the ground 
state of Ne to the Rydberg series ls^2s2p^np with n = 3 to 8. 



Resonant transition 


Approach 


CTO {Mb) 


r (meV) 


Er (eV) 


q 


2 


2s 3p 


TDSLHF" 


8.34 


12.80 


43.358 


-4.03 


0.551 




S-TDSLHF'' 


8.40 


15.23 


45.447 


-3.40 


0.619 




S-EXX/ALDA'^ 


8.09 


17.9 


45.438 


-3.18 


0.547 




S-TDLSDA'^ 


8.28 


18.04 


45.453 


-2.40 


0.764 




TDLDA'^ 


8.18 


13.90 


46.253 


-3.69 


0.514 




R-matrix^ 




34.9 


45.534 








Expt.9 


8.6 ±0.6 


13 ± 2 


45.546 ± 0.008 


-1.6 ±0.2 


0.70 ±0.07 


2s Ap 


TDSLHF" 


8.12 


3.89 


45.028 


-4.35 


0.543 




S-TDSLHF'' 


8.15 


4.55 


47.115 


-3.74 


0.606 




S-EXX/ALDA'^ 


7.89 


5.5 


47.093 


-3.35 


0.528 




S-TDLSDA'* 


8.06 


5.14 


47.098 


-2.62 


0.783 




TDLDA<= 


7.98 


3.86 


47.397 


-3.95 


0.505 




R-matrix-^ 




6.65 


47.111 








Expt.9 


8.0 ± 0.6 


4.5 ± 1.5 


47.121 ± 0.005 


-1.6 ± 0.3 


0.70 ± 0.07 


2s 5p 


TDSLHF" 


8.05 


1.55 


45.625 


-4.47 


0.537 




S-TDSLHF'' 


8.06 


1.81 


47.713 


-3.85 


0.600 




S-TDLSDA'* 


7.91 


2.20 


47.683 


-2.72 


0.783 




TDLDA'^ 


7.91 


1.62 


47.814 


-4.05 


0.502 




R-matrix^ 




2.47 


47.692 








Expt.f 


8.2 ±0.6 


2±1 


47.692 ± 0.005 


-1.6 ±0.5 


0.70 ±0.14 


2s — > 6p 


TDSLHF" 


8.00 


0.89 


45.901 


-4.57 


0.539 




S-TDSLHF'' 


8.00 


1.04 


47.989 


-3.94 


0.601 




R-matrix^ 




1.28 


47.967 








Expt.9 






47.967 ± 0.006 






2s 7p 


TDSLHF" 


7.99 


0.43 


46.075 


-4.99 


0.555 




S-TDSLHF'' 


8.01 


0.59 


48.163 


-3.88 


0.585 




R-matrix^ 




0.73 


48.119 








Expt.9 






48.116 ±0.006 






2s 8p 


TDSLHF" 


7.94 


0.37 


46.171 


-4.54 


0.534 




S-TDSLHF'' 


7.96 


0.42 


48.258 


-3.92 


0.605 




R-matrix^ 




0.46 


48.211 








Expt.f 






48.207 ± 0.006 







" TDSLHF results. TDSLHF results with the experimental 2s electron orbital energy. " EXX/ALDA results witji the 
experimental 2s electron orbital energy jlQ|. TDLSDA results with the experimental 2s electron orbital energy 
TDLDA results j]. ^ R- matrix results ^ experimental results [s^]. 
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TABLE II: Autoionization resonance peak positions (in eV) for the photoionizations from Ne outer-shell 
excited states < ls|ls^2s|2s|2p?2p? > n/f with n = 3 4 and Z = 2. 



Resonant 
transition 




Initial 






states 




{}*3st 


{}3PT 


{}3dt 


{}4ST 


{}4PT 


{}4dT 


2pi- — »3si- 


- 


23.829 


26.150 


23.527 


24.863 


26.509 


4si- 


24.526 


26.738 


29.139 




27.693 


29.554 


5sf 


25.383 


27.630 


30.153 


26.839 


28.581 


30.452 


6sf 


25.759 


28.014 


30.575 


27.206 


28.972 


30.896 


7st 


26.013 


28.221 


30.803 


27.405 


29.175 


31.119 


8sf 




28.406 


30.928 


27.549 


29.293 


31.250 


9sf 






31.043 




29.386 


31.334 


lOsf 












31.416 


ad^ 


24.855 


27.078 


29.565 


26.371 


28.060 


30.001 


4df 


25.516 


27.769 


30.330 


26.999 


28.730 


30.670 


5d| 


25.826 


28.085 


30.670 


27.290 


29.046 


30.986 


6d| 


26.066 


28.264 


30.850 


27.453 


29.217 


31.173 


7d| 






30.961 


27.590 


29.320 


31.282 


8d| 










29.410 


31.358 


2pi -+3s| 


29.984 


30.382 


30.656 


30.594 


30.587 


30.643 


4s| 


35.446 


36.948 


38.389 


38.224 


38.302 


38.559 


5s| 


36.809 


38.447 


40.406 


40.463 


40.722 


41.081 


6s| 










42.194 


42.635 


3d; 


36.036 


37.378 


38.327 


38.161 


38.213 


38.412 


4d; 


37.018 


38.580 


40.479 


40.406 


40.697 


41.021 


5d; 








41.693 


41.995 


42.423 


2sf — +2pi- 


22.994 


23.032 


23.032 


23.048 


23.081 


23.056 




46.964 


49.187 


51.685 


48.719 


50.243 


52.129 


4p| 


48.510 


50.798 


53.310 


50.044 


51.767 


53.682 


5PT 


49.075 


51.381 


53.960 


50.572 


52.352 


54.265 


6PT 


49.359 


51.659 


54.260 


50.833 


52.630 


54.577 


'PT 


49.623 


51.835 


54.420 


50 986 


52.782 


54.745 


8PT 




52.036 


54.532 


51.136 


52.883 


54.845 


9PT 






54.654 




52.978 


54.923 


lOpT 










53.098 


55.013 


llPT 












55.114 


2s| ^3p| 


56.790 


57.457 


57.892 


57.786 


57.797 


57.911 


4p| 


59.650 


61.242 


62.926 


62.818 


62.940 


63.218 


5p; 


60.817 


62.458 


64.412 


64.404 


64.747 


65.147 



* {} = ls|ls|2s|2sf 2pj2p^ is the abbreviation of the core electron configuration. 
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TABLE III: Autoionization resonance peak positions (in eV) for the photoionizations from Ne inner-shell 
excited states < ls|lsx2s|2p?2p? > nli with n = 3 ~ 4 and Z = ^ 2. 



Resonant 
transition 




Initial 






states 




{}*3si 


{}3Pi 


{}3di 




{Hpi 


{}4d| 


2p| -»3s| 




21.655 


23.821 


21.478 


22.681 


24.128 


4s 


22.188 


24.648 


26.738 




25.356 


27.086 


5s 1 


23.026 


25.595 


27.742 


24.626 


26.245 


27.973 


6s 1 


23.394 


25.998 


28.167 


24.994 


26.656 


28.403 


7si 


23.644 


26.215 


28.381 


25.187 


26.866 


28.629 


8s| 




26.411 


28.509 


25.328 


26.988 


28.756 


9si 






28.624 


25.499 


27.094 


28.841 


lOs; 






28.768 




27.224 


28.947 


3d| 


22.488 


24.923 


27.141 


24.163 


25.739 


27.508 


4d; 


23.143 


25.720 


27.908 


24.784 


26.398 


28.150 


5di 


23.453 


26.060 


28.248 


25.073 


26.730 


28.493 


6di 


23.686 


26.254 


28.427 


25.233 


26.907 


28.678 


7di 




26.444 


28.542 


25.369 


27.015 


28.790 


8di 






28.654 


25.543 


27.118 


28.866 


9di 






28.801 




27.249 


29.048 


2p| — >3s| 


30.711 


30.942 


31.198 


31.122 


31.135 


31.201 


4s I 


36.393 


37.647 


39.089 


38.904 


39.002 


39.263 


5s f 


37.780 


39.152 


41.141 


41.176 


41.440 


41.821 


6s f 










42.937 


43.440 


3d| 


36.994 


38.088 


39.037 


38.842 


38.918 


38.108 


4dT 


37.987 


39.288 


41.206 


41.127 


41.421 


41.750 


5d| 






42.412 


42.390 


42.725 


43.190 


2si- ^3p| 


57.854 


58.344 


58.768 


58.638 


58.676 


58.785 


4p| 


60.815 


62.154 


63.862 


63.729 


63.865 


64.151 


5PT 


61.998 


63.375 


65.351 


65.337 


65.686 


66.107 



* {} = ls|ls|2s|2p^2p| is the abbreviation of the core electron configuration. 
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